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Absolute  instability  of  a  liquid  jet  in  a  gas 

S.  P.  Lin  and  Z.  W.  Lian 

Department  of  Mechanical  and  Industrial  Engineering.  Clarkson  University,  Potsdam.  Sew  York  13676 

(Received  15  September  1988;  accepted  8  November  1988) 

effect  of  the  ambient  gas  density  on  the  onset  of  absolute  instability  in  a  viscous  liquid  jet  is 
examined.  The  critical  Weber  number,  above  which  the  instability  is  convective  and  below 
which  the  instability  is  absolute,  is  determined  as  a  function  of  Reynolds  number  and  the 
density  ratio  of  gas  to  liquid.  It  is  shown  that  the  gas  density  has  the  effect  of  raising  the 
critical  Weber  number.  It  also  raises  the  cutoff  wavenumber  below  which  disturbances  are 
spatially  amplified  and  above  which  they  are  damped.  |  ' 


I.  INTRODUCTION 

"  The  capillary  instability  of  an  infinitely  long  jet  with 
respect  to  temporally  growing  disturbances  was  analyzed  by 
Rayleigh?^  Keller’  etal.  examined  thecapillary  instability  of 
a  semi-infinite  jet  with  respect  to  spatially  growing  distur¬ 
bances.  They  found  that  the  temporal  and  spatial  distur¬ 
bances  are  analytically  related  if  the  Weber  number  is  suffi¬ 
ciently  large.  For  sufficiently  small  Weber  numbers,  Leib 
and  Goldstein^ound  that  the  state  of  convective  instability 
obtained  by  Keller  et  al.  actually  cannot  be  reached  by  a 
given  initial  disturbance  in  the  sense  of  Briggs’^  and  Bers  * 
Leib  and  Goldstein  demonstrated,  for  the  first  time,  the  exis¬ 
tence  of  the  absolute  instability  in  an  inviscid  jet.  Recently, 
they'’  determined  from  Chandrasekhar’s'  dispersion  equa¬ 
tion  the  critical  Weber  number  below  which  a  viscous  jet  is 
absolutely  unstable  as  a  function  of  Reynolds  number.  They 
also  found  that  the  cutoff  wavenumber,  above  which  the  dis¬ 
turbance  is  spatially  damped,  is  one  independent  of  the  We¬ 
ber  number. 

Here,  we  examine  the  effect  of  the  ambient  gas  density 
on  the  absolute  instability  discovered  by  Leib  and  Goldstein. 
It  is  shown  that  the  gas  density  has  the  effect  of  enlarging  the 
domain  of  absolute  instability  in  the  Reynolds- Weber  num¬ 
ber  space.  Moreover,  it  raises  the  cutoff  wavenumber  below 
which  the  jet  is  convectively  unstable,  and  also  raises  the 
amplification  rate  of  the  spatially  growing  disturbances 
when  the  absolute  instability  is  absent. 

II.  FORMULATION 

Consider  a  circular  cylindrical  jet  of  an  incompressible 
viscous  Newtonian  liquid  issued  from  a  nozzle  into  an  un¬ 
bounded  inviscid  gas.  The  governing  dynamic  equations  of 
motion  in  the  liquid  and  the  gas  phases  are,  respectively,  the 
Navier-Stokes  and  Euler  equations.  The  boundary  condi¬ 
tions  are  the  vanishing  of  the  net  force  per  unit  area  of  the 
interface,  and  the  equality  of  radial  fluid  velocity  in  each 
phase  with  the  total  time  rate  of  change  of  the  interfacial 
position.  A  uniform  velocity  distribution  f/  in  a  circular  jet 
of  radius  r„  in  a  quiescent  gas  is  an  exact  solution  to  this  set  of 
differential  equations  in  the  absence  of  gravity.  This  exact 
mathematical  solution  representing  a  possible  basic  state  is 
physically  unstable.  The  stability  analysis  of  this  basic  state 
with  respect  to  any  Fourier  component  of  the  disturbances 


of  the  form  Q  exp(iur  -I-  iky)  led  to  the  characteristic  equa¬ 
tion  ( 1 )  when  viewed  in  a  reference  frame  with  its  origin 
fixed  at  the  nozzle  exit,"  ’ 


(<u  -  ik)^  + 


2k^  n\ 

R\lo 


(k) 


Xlct)  —  ik)  +  CO 


o(^) 

2P2 


2kA  I,{k)I\{A)^ 
+  k'^  Io(k)I,U)i 
-k^  Ko(k)f(k) 


Pi  +  k'^  K,(k)Io(k) 
k^  I,{k) 


_-L^(l  _yt2)  i4 

/?  A‘  +  k^  J„(k) 


=  0, 


(1) 


where  C„  is  the  wave  amplitude,  co  is  the  complex  wave  fre¬ 
quency  { the  real  part  of  which  gives  the  exponential  tempo¬ 
ral  growth  rate,  and  imaginary  part  of  which  gives  the  wave 
frequency  of  disturbances),  k  is  the  complex  wavenumber 
(the  real  part  of  which  is  equal  to  2trr,/(  wavelength)  and 
the  imaginary  part  of  which  gives  the  spatial  amplification 
rate),  r  is  time  normalized  with  r^/U,  y  is  the  axial  distance 
measured  in  the  unit  of  r„  in  the  opposite  direction  of  the  jet 
flow,  /  and  K  are,  respectively,  the  modified  Bessel  func¬ 
tions  of  the  first  and  second  kind,  their  subscripts  denote  the 
order  of  the  functions,  and  A  is  defined  by 


A'  =  k^  +  Rico  —  ik). 


The  three  independent  flow  parameters  in  ( 1 )  are  Reynolds 
number  R,  Weber  number^,  and  the  density  ratio  Q  defined 
respectively,  by 

R=Ur^v„  0  =  pyU^ra/a,  Q  =  P2/p„ 

where  v  is  kinematic  viscosity,  p  is  density,  a  is  surface  ten¬ 
sion,  and  the  subscripts  1  and  2  denote,  respectively,  the 
liquid  and  gas  phases.  Note  that  the  above  nondimensionali- 
zation  of  time  implies  that  the  dimensional  frequency  is  giv¬ 
en  by  wf/  /rg. 


III.  RESULTS 

Figure  I  shows  the  loci  of  the  two  characteristic  roots 
obtained  when  a  given  set  of  flow  parameters  as  the  temporal 
growth  rate  co,  varies  for  a  number  of  values  of  the  wave 
frequency  co,.  As  the  co,  is  reduced  from  a  finite  positive 
value  to  zero,  these  loci  approach  toward  two  different 
branches  of  amplification  curves  co,  =0  for  the  spatially 
growing  disturbances.  Note  that  k,>0  gives  the  spatial 
instability  because  the  y  axis  is  chosen  to  be  negative  in  the 
flow  direction.  It  can  be  seen  that  two  of  the  loci  with 
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0.74  <<u,  <0.75  will  intersect  in  the  upper  half-plane  at  a  gives  the  dependence  of  the  critical  Weber  number  on for 

saddle  point.’  A  similar  situation  has  already  been  shown  by  three  values  of  Q.  Here,  Q  =  0.00 1 3  corresponds  to  the  air  to 

Leib  and  Goldstein’  for  the  case  of  negligible  gas  density.  As  water  density  ratio  in  one  atmospheric  pressure  in  room  tem- 

the  Weber  number  is  increased,  with  the  values  of  Q  and  It  perature.  The  curve  Q  =  0  was  obtained  by  Leib  and  Gold- 

given  in  Fig.  1  being  fixed,  this  saddle  point  moves  closer  to  stein.  In  the  parameter  range  below  the  curves,  the  jet  is 

the  pure  spatial  amplification  curve  for  disturbances  with  a  absolutely  unstable.  The  dashed  line  indicates  the  inviscid 

“group  velocity”  in  the  direction  of  jet  flow.  Note  that  the  limit  obtained  by  Leib  and  Goldstein  for  the  case  of  0  =  0.  It 

amplification  curve,  which  does  not  originate  from  k,  =  0,  is  is  seen  that  the  gas  density  has  the  effect  of  enlarging  the 

for  the  disturbances  with  upstream  propagating  “group  ve-  domain  of  absolute  instability  for  the  0-R  plane, 

locity.”  Following  Leib  and  Goldstein,  the  minimum  value  Above  these  curves,  the  jet  is  convectively  unstable.  Fig- 

of  P  for  which  the  saddle  point  of  Eq.  ( 1 )  lies  in  the  upper  ures  3  and  4  give  a  few  typical  amplification  curves  for  con- 
half-plane  is  defined  as  the  critical  Weber  number.  Figure  2  vectively  unstable  disturbances.  Leib  and  Goldstein  showed 


FIG.  3.  The  growth  rate  for  convec- 
tively  unstable  disturbances: 
/!  =  20  000,^=400. 


that  the  cutoff  wavenumber,  above  which  the  disturbance  is 
damped,  is  one  independent  of  P  for  the  case  of  Q  =  0.  It  is 
seen  from  Figs.  3  and  4  that  the  gas  density  tends  to  raise  the 
cutoff  wavenumber  very  significantly.  This  has  a  practical 
significance  of  producing  smaller  droplets  by  use  of  larger 
ambient  gas  density.’ 

It  is  very  hard  to  resist  the  temptation  of  speculating 
that  the  absolute  instability  may  correspond  to  the  dripping 
mode  of  the  jet  instability.  When  a  saddle  point  of  the  eigen¬ 
value  exists  at  a  certain  frequency  and  there  exist  both  up¬ 
stream  and  downstream  propagating  convectively  unstable 
branches,  a  disturbance  with  this  certain  frequency  may  be 
alternatively  attracted  toward  these  two  branches  as  time 


evolves.  The  simultaneous  spatial  and  temporal  growth  al¬ 
ternating  in  the  upstream  and  downstream  directions  with  a 
regular  frequency  casually  observed  in  a  dripping  jet  seems 
to  correspond  to  the  mathematical  picture  of  absolute  insta¬ 
bility.  This  speculation  is  even  more  tempting  when  one  real¬ 
izes  that  the  absolute  instability  is  predicted  only  for  a  small 
Weber  number,  which  is  the  ratio  of  the  inertial  force  to  the 
surface  tension  force.  However,  the  argument  remains  spec¬ 
ulative,  since  the  dripping  phenomenon  involves  a  highly 
nonlinear  finite  amplitude  of  disturbances  that  cannot  be 
accounted  for  in  the  present  linear  theory.  However,  Leib 
and  Goldstein^  suggested  that  the  absolutely  unstable  distur¬ 
bance  may  grow  nonlinearly  by  the  mechanism  of  direct  res- 
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onance  advanced  by  Akylas  and  Benney. The  physical  sig¬ 
nificance  of  the  absolute  instability  in  the  context  of 
nonlinear  theories  remains  to  be  explored. 
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